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Abstract 

Many signals in Nature, technology and experiment have a multi- 
5j ■ component structure. By spectral decomposition and projection on 

5a . the eigenvectors of a family of unitary operators, a robust method is 

developed to decompose a signals in its components. Different signal 
traits may be emphasized by different choices of the unitary family. 
The method is illustrated in simulated data and on data obtained from 
plasma reflectometry experiments in the tore Supra. 



> 

Jj : 1 Introduction 

(^! , Most natural and man-made signals are nonstationary and may be thought 

of as having a multicomponent structure. Bat echolocation, whale sounds, 
radar, sonar and many others are examples of this kind of signals. The notion 
of nonstationarity is easy to define. However, the concept of signal component 
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is not so clearly defined. Because time and frequency descriptions are stan- 
dard methods of signal analysis, many authors have attempted to base the 
characterization of signal components on the analysis of the time-frequency 
plane. There is a large class of time- frequency signal representations (TFR). 
An important set of such TFR's is Cohen's classp], obtained by convolutions 
with the Wigner distribution 

C $ (*,/) = / fw {u, v)${t-uj-v) dudv (1) 

W (u, v) being the Wigner distribution 

W(tJ) = Js(t + I) s* (t - I) e-^f dT (2 ) 

Once one particular TFR of the signal is constructed, the search for com- 
ponents may be done by looking for amplitude concentrations in the time- 
frequency plane. This is the methodology that has been followed by most 
authors^] [31 0] [S] [H] H [H] IH] HO] [H] - The notions of instantaneous fre- 
quency and instantaneous bandwidth play an important role in these studies. 

An important drawback from the use of TFR's is the fact that they may 
have negative terms, cross terms or lack the correct marginal properties in 
time and frequency. Even if, by the choice of a clever kernel or a smoothing 
or filtering operation, the TFR's are apparently free from these problems, 
there is no guarantee that they are free from artifacts that might lead to 
unwarranted inferences about the signal properties. This is a consequence 
of the basic fact that time (t) and frequency (u> = zjj), being associated to 
a pair of noncommuting operators, there can never be a joint probability 
distribution in the time-frequency plane. 

Our approach to component separation starts from the insight that the 
notion of component depends as much on the observer as on the observed 
object. That is, when we speak about a component of a signal we are in 
fact referring to a particular feature of the signal that we want to emphasize. 
For example, if time and frequency are the features that interest us, it might 
indeed be the salient features in the time-frequency plane to be identified as 
components. However, if it is frequency and fractality (scale) that interests 
us, the notion of component and the nature of the decomposition would be 
completely different. 

In general, the features that interest us correspond to incompatible no- 
tions (that is, to noncommuting operators). Therefore to look for robust 



characterizations in a joint feature plane is an hopeless task because the 
noncommutativity of the operators precludes the existence of joint probabil- 
ities. Instead, in our approach, we consider spectral decompositions using 
the eigenvectors of linear combinations of the operators. The sum of the 
squares of the signal projections on these eigenvectors has the same norm as 
the signal, thus providing an exact probabilistic interpretation. Important 
operator linear combinations are the time-frequency 



B (s) (jj,, v) = jit + vuj = jit + vi— (3) 



the frequency-scale 



B[ A) (/2, v) = jjuj + vD = fMV + v- (tto + tot) (4) 

and the time-scale 

B { 2 A) (/!, u)=fit + vD (5) 

Then, a quadratic positive signal transform is defined by 

Mf (X, n,v)= f s* (t) 5 (B (n, v)-X)s (t) dt (6) 

called a B— tomogram which, for a normalized signal 

\s(t)\ 2 dt=l (7) 

is also normalized 

/ Mf (X, /i, v)dX = l (8) 

For each (/i, v) pair, the tomograms Mf (X, //, v) provide a probability dis- 
tribution on the variable X, corresponding to a linear combination of the 
chosen operators (time and frequency, frequency and scale or time and scale). 
Therefore, by exploring the family of operators for all pairs (/x, v) one obtains 
a robust (probability) description of the signal at all intermediate operator 
combinations. 

Using the (symmetric) operators B (/i, v) and their corresponding unitary 
exponentiations 

U (/j,, v) — exp(iB (/j,, v)) (9) 



a unified description of all currently known integral transforms has been 
obtained [1 2] . Explicit expressions for the tomograms in the three cases flMHD 
may be found in |13j . 

Of particular interest for the component analysis in this paper is the 
time-frequency operator B^ (/i, v) for which 

2 



M s (x, fjL, V) 



I s(t)exp ( — t 2 - — t J dt 



(10) 



2ir\v\ 
called the symplectic tomogram. The tomogram is an homogeneous function 

M 3 (-A-) =\p\M.(x,»,u) (11) 



p p p ^ 

For the particular case \x = cos9,v = sin#, the symplectic tomogram 
coincides with the Radon transform [14j. already used for signal analysis by 
several authors [15] [16] [T7] in a different context. 

Once a tomogram for B = fiO\ + Z/O2 is constructed, what one obtains 
in the (X, (/i, v)) (hyper-) plane is an image of the probability flow from the 
0\— description of the signal to the O2— description, through all the interme- 
diate steps of the linear combination. In contrast with the time-frequency 
representations we need not worry about cross-terms or artifacts, because of 
the exact probability interpretation of the tomogram. Then, we may define 
as a component of the signal any distinct feature (ridge, peak, etc.) of the 
probability distribution in the (X, (/i, v)) (hyper-) plane. It is clear that the 
notion of component is contingent on the choice of the pair (0 1; 2 )- 

In Sect. 2 we analyze in detail the time-frequency tomogram, the choice 
of a complete orthogonal basis of eigenvectors of B s (/1, v) for the projection 
of the signal and how the component identification may be carried out by 
spectral decomposition into subsets of this basis. In Sect. 3 a few examples 
of component decomposition of noisy signals are worked out, which show 
the effectiveness of the method. Finally, in Sect. 4 the method is applied to 
experimental data obtained in the reflectometry analysis of plasma density. 
In the Appendixes we collect a few results, which are useful for the practical 
calculation of the symplectic tomograms. 

2 Tomograms and Signal Analysis 

Following the ideas described in the introduction, a probability family of 
distributions, M s (x,9), is defined from a (general) complex signal s(t), t G 



[0,T]by 



with 



M a { x ,e) 



s(t)* e /(t)dt 
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<s^/> (12) 



TflT/ . 1 /—icosO 9 ix \ , 

q e /(t) = —exp( , - t 2 + — —t) 13 

y/T \ 2 sin 9 sin 9 ) 

This is a particular case of Eg. ([TO]) for \x = cos9,u = sin9. Here 9 is a 
parameter that interpolates between the time and the frequency operators, 
thus running from to n/2 whereas x is allowed to be any real number. Notice 
that the ^' T 's are generalized eigenfunctions for any spectral value x of the 
operator U(9). Therefore M s (x, 9) is a (positive) probability distribution as a 
function of x for each 9. From an abstract point of view, since for different 0's 
the U{9) (see Eq.(j9])) are unitarily equivalent operators, all the tomograms 
share the same information. However, from a practical point of view the 
situation is somehow different. In fact when 9 changes from to tt/2 the 
information on the time localization of the signal will gradually concentrate 
on large x values which are unattainable because of sampling limitations. 
On the other hand and by opposite reasons, close to 9 = the frequency 
information is lost. Therefore we search for intermediate values of 9 where 
a good compromise may be found. For such intermediate values, as we shall 
see in several examples, it is possible to pull apart different components of 
the signal that take into account both time and frequency information. The 
reason why this is the case will be clear by looking at the properties of (TLB"]) . 
First we select a subset x n in such a way that the corresponding family 
|\l/^'J(t)} is orthogonal and normalized, 

< \V e ' T \V e ' T >= 5 (14) 

This is possible by taking the sequence 

x n = x + -=r sin 9 (15) 

where Xq is freely chosen (in general we take Xq = but it is possible to 
make other choices, depending on what is more suitable for the signal under 
study). 

A glance at the shape of the functions fflBl shows that the nodes (the 
zero crossings) t n of the real (resp. imaginary) part of ^ e x ,T are the solutions 



of 

cos x 

ir—jl-—Jn = 1™ {resp.27in + 7i/2) (16) 

2 sin 9 sin 9 

It is clear that |t n +i — t n \ scales as y/n and that, for fixed 9, the oscillation 
length at a given t decreases when \x\ increases. As a result, the projection 
of the signal on the |\l/^(t)} will locally explore different scales. On the 
other hand, changing 9 will modify the first term of (TTBT) in such a way that 
the local time scale is larger when 9 becomes larger in agreement with the 
uncertainty principle. 

We then consider the projections of the signal s(t) 

cl(s)=<s,^> (17) 

which in the following are used for signal processing purposes. In particular 
a natural choice for denoising consists in eliminating the c e x (s) such that 

KM 2 <e (18) 

for some chosen threshold e, the remainder being used to reconstruct a de- 
noised signal. In this case a proper choice of 9 is an important issue in the 
method. 

In the present work we mainly explore the spectral decomposition of the 
signal to perform a mult i- component analysis. This is done by selecting 
subsets Tk of the x n and reconstructing partial signals (/c-components) by 
restricting the sum to 

«*(*) = E 4w*ST (*) ( 19 ) 

for each k. 

Eq. (fl9|) builds the signal components as spectral projections of s. As we 
shall see, by an appropriate choice of 9, it is possible to use this technique 
to disentangle the different components of a signal. This work is mainly 
devoted to the analysis of reflectometry signals in plasma physics for which 
the method seems well adapted. 

Notice that the inverse of \t e n+l — t e n \ plays the role of a quasi- instantaneous 
frequency defined in a #-like scale. This is a piecewise constant function but, 
as seen from ffl3|) . it grows approximately linearly in time with slope tan _1 (0). 
We used such time scales to control the quality of the sampling. 

6 



3 Examples: Simulated data 

In this section we discuss the general method presented in the previous sec- 
tion in two particular simulated signals. The first example shows how the 
method is able to disentangle a signal with different time as well as frequency 
components. In the second example a signal with time-varying frequency is 
analyzed. 

3.1 First Example 

Let us consider a signal y(t), of duration T = 20 s, that is the sum of three 
sinusoidal complex signals y^, k = 1,2, 3, plus a noise component b : 

y(t)=y 1 (t)+y 2 (t)+y 3 (t)+b(t) (20) 

where 

yi it) = exp (i25t) , t G [0, 20] 
y 2 (t) = exp (i75t) ,te[0, 5] 
y 3 (t) = exp{i75t),t G [10,20] 

The Signal to Noise Ratio, SNR y ^, is about 10 dB, the SNR being defined 
by 

SNR(y,b) = 101og 10 § (21) 

with P y = Tf, J \y(t)\ dt and P\, = ^ f Q \b(t)\ dt. The real part of the 
simulated data, 1Z [y(i)], is shown in the FigJH In order to test the robustness 
of the projection protocol we first compare the original signal y(t) with a 
reconstructed signal y(t) given by 

175 

m= E 4d/)*Sf(*) ( 22 ) 

x n =-175 

The quadratic error E(y, y), between the original and the reconstructed signal 
is less than —27dB. The quadratic error is defined as : 

£(?/,£) = 10 log 10 %^ (23) 




Figure 1: Temporal representation of the simulated data TZ [y(t)) defined by 
equation fl20|) . 



The quadratic error grows to —22dB if the reconstruction is limited to the 
range used for the component analysis below (i.e. 45 < x n < 155). 

The value 9 = f is chosen by direct inspection of the tomogram of the 
signal y(t). In this case one sees three well differentiated spectral components 
(Figj2]). Clearly this is not the unique possible (even logical) choice. From 
a practical point of view, either some preliminary information is available 
on the structure of the signal, or we may try different choices of 9 knowing 
that the incertitude on the time support will increase with 9 whereas the 
quasi-local frequency incertitude will decrease. 

In this example, we performed the factorization of y(t) in three compo- 
nents y~i(t), ij2{t)i and yzif) defined respectively by the equations (1241 . (1251) 
and (1261) . Using different values of 9, the quadratic errors E(y, y), E(yi, y~i), 
E(y 2 ,y2) and E(y 3 ,y 3 ) are computed (Eq. [23]). We summarize the corre- 
sponding data in the following table: 



9 


vr/8 


tt/5 


3tt/10 


4vr/5 


tt/2 


E(yi,m) 


-14.5dB 


-17.5dB 


-18.5dB 


-17.5dB 


-12.5dB 


£(2/2,2/2) 


-10.5dB 


-12.5dB 


-9dB 


-7dB 


-0.5dB 


E(y 3 ,y 3 ) 


-14.5dB 


-14dB 


-13.5dB 


-7dB 


-4dB 


E(y,y) 


-26.5dB 


-27dB 


-30dB 


-30dB 


-28dB 



In cases where no a priori information is available, different choices of 



9 may provide meaningful information about the signal structure. In this 
example, by looking at the data presented on the table above, the choice 
of 9 — | to carry out the factorization seems to give the best performance. 
Then we simply apply an energy threshold e = 0.1, which is about 15% of 
the energy level of the signal, to decompose the signal y in three components 
(Fig©. 
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Figure 2: c e Xn {y) spectrum of the simulated data y(t) for 9 = | and < x n < 

175 



The first component, y~\{t) corresponds to the spectral range 45 < x n < 



65: 



65 



m= E oxrw 



(24) 



£n=45 



The second component, 2/2 (£) corresponds to the spectral range 135 < 

115 



x n < 145: 



x n =135 



(25) 



The real part of 1/2 (t) and y 2 (t) are presented FigJ3 

The last component, y^{t) corresponds to the spectral range 145 < x n < 



155: 



155 



y?> 



(*)= E c6 *MWAt) 



(26) 



x„=145 



FigS] gives a representation of both TZ [2/3 (t)] and TZ [2/3 (t)]. 
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time 



Figure 3: 1Z[y 2 (t)] and TZ[y 2 (t)]. For visual purposes, the mean value of 
1Z [2/2^)] is shifted to 3. 

The quadratic errors E(yi,y~i), E(y2,y~2) and E(y 3 ,y 3 ) can be read from 
the table above. They are, respectively, — 17.5dB, —12.5dB and — lAdB. 

For comparison, the projection of the simulated data y(t) in the frequency 
domain (#0 = f ), presented on figure [5j show that the factorization in three 
components is not possible : only two components can be extract from this 
projection. At the frequency x n = 25rd/s, the component will be equal to 
y~i(t). At the frequency x n = 75rd/s, it will be impossible to set apart y2(f) 
and 2/ 3 (t) and the component will be equal to y 2 it) + y~zit). 

3.2 Second Example 

Here we analyze the decomposition into elementary components of another 
signal which aims to mimic, in a simplified way, the case of an incident plus a 
reflected wave delayed in time and with an acquired time-dependent change 
in phase. In this case the simulated signal y(t) is the sum of an "incident" 
chirp yo(t) and a "deformed reflected" chirp ynit). White noise is added to 
the signal. The incident chirp is: 



yo(t) 



a »*o(t) 



(27) 



with $ (t) = a t 2 + b t. 

The "instantaneous frequency" of yoit) sweeps linearly from 75 rd/s to 50 
rd/s during 20s. Its phase derivative is linearly dependent on time :^$o(^0 = 
2a t + b . 
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Figure 4: Representation of 7?. [2/3 (t)] and 1Z [y~3(t)]. For visual purposes, the 
mean value oflZ [2/3 (t)] is shifted to 3. 

The "reflected" signal ?/_r(£) is delayed by £# = 3s from the incident one 
and continuously sweeps from 75rd/s to 50r<i/s : 



2/fl(*) 



,**«(*) 



(28) 



where $#(£) = a#(£ — £r) 2 + 6r(£ — £r) + 10(t — £r) 5 - In this case the phase 
derivative ^$^(t) is not a linear function. This signal is zero during the first 
3s seconds and ends up at t = 23s. 
The simulated signal is defined by : 



y{t)=yo{t) + y R {t) + b{t) 



(29) 



The Signal to Noise Ratio, SNR(y,b), is equal to lbdB. The real 1Z[y(t)) 
part of this signal is shown in Fig^J 

FigfTJshows 4<&o(£) and 4$i?(t) as a function of time. Notice that, except 
for the three first seconds, the spectrum of the signals yo(t) and yuif) is almost 
the same. 

The tomogram of the first 20s of y(t), M y (8,x) = < y, ^^ > has 
a maximum for sin(0) « 0.6 (FigJHJ) corresponding to the "incident" part 
of the signal that mainly projects in the unique \E^' T that matches $o(^)- 
Then, we take the value of 9 — n/5 to carry out the separation of y(t) in its 
components. 

The corresponding spectrum c e Xn (y) is shown in FigJH Based on this 
spectrum we decompose the signal in two spectral components. 
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Figure 5: Projection c e x ° of the signal for 9q 



From the first component we reconstruct the "incident" chirp yo(t) by: 



47.25 



Ut)= £c*M*2f(*) 



(30) 



x n =45 



The quadratic error, between y~o(t) and yo(t), E(y ,y ), is —9.5dB. 

From the second spectral component we reconstruct the "reflected" chirp 
given by: 

50.5 

y R (t)= J2 i»**J(0 ( 31 ) 

x„=47.5 

In this case the quadratic error E(y R ,y R ) is — lOdB. This may be compared 
with a quadratic error E(y, y) of —29dB for the total signal reconstructed 
from the spectral projection corresponding to 45 < x n < 50.5. 

We have tested the method with different delays that encode for the 
distance of the two "frequencies" . The quality of the disentanglement deteri- 
orates when the delay decreases. But it can still be done for a delay as short 
as Is. Therefore we conclude that the method is quite robust. 
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Figure 6: Temporal representation of the simulated data TZ [y(t)) defined by 
equation fl29|) . 



4 An application to reflectometry data. Com- 
ponent analysis 

We now use signals coming from reflectometry measurement in plasma physics, 
to show the ability of the tomogram methods to separate different compo- 
nents of the signal to which it is then possible to assign a clear physical mean- 
ing. The reflectometry diagnostic is widely used to determine the electronic 
density profile in a tokamak. The principle, based upon a radar technique 
|18j . is to measure the phase of a probing wave reflected by the plasma cut-off 
layer at a given density, where the refractive index goes to zero. The deter- 
mination of the density profile can be achieved by continuously sweeping the 
frequency of the probing wave. 

Different techniques are used to measure the density profile on fusion 
plasmas [19] (phase difference, ultrashort pulses, continuous sweep, ...). A 
broadband reflectometer operating in the frequency range 50-75 GHz (V 
band) [20], [21] and 75-110 GHz (W band) [22] has been developed on Tore 
Supra to measure the electron density profiles at the edge. 

The sweep frequency reflectometry system of Tore Supra launches a prob- 
ing wave on the extraordinary mode polarization (X mode) in the V band 
(50-75 GHz). The emitting and receiving antennas are located at about 1.20 
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Figure 7: Representation of ^$ (^) and Jj$#(£) as a function of time. 

m from the plasma edge outside the vacuum vessel. The reflectometry sys- 
tem operates in burst mode, i.e. the sweeps are performed repeatedly every 
25/is. The duration of one sweep, E (t) = Aoe % ^\ is 20/is and 5000 chirps 
are send during one measurement. During the 20/xs measurement time, the 
frequency of the probing wave is continuously sweeping, from 50 GHz to 75 
GHz (V band). 

The heterodyne reflectometers, with I/Q detection, provide a good Signal 
to Noise Ratio, up to AOdB. For each sweep, the reflected chirp E R (t) is mixed 
with the incident sweep Eo(t) and only the interference term is recorded as 
in-phase and 90° phase shifted signals sampled at T e = 10~ 8 s 

xi(t) = A A R (t) cos(<p(t)) 

x 2 (t) = A A R (t) sm(ip(t)) 
For each sweep, the phase <p(t) of the reflected signal is represented by 

y{t) = Xl (t) + ix 2 (t) = A{t)e iip{t) (32) 

The amplitude of this signal A(t) = A A R (t) is low frequency. The real part 
of one such signal y(t) is shown in FigJTOl 

The tomogram M y (x, 9) of the signal is shown in FigfTTlwhere it is possible 
to see that it carries three main different components. The choice of 9 = |^ 
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Figure 8: Tomogram of the simulated data defined by Eq. fl29l) 



to perform the factorization of the signal was done by simple inspection of 
this tomogram. The spectrum c e Xn (y) of the reflectometry signal for 6 = || 
is shown in Figfl2l 

When reconstructing y(t) by: 



200 



m= E 4(»)^) 



(33) 



=-200 



the quadratic error E(y,y), between the original and the reconstructed sig- 
nals is —25dB. 

Factorization of the reflectometry signal 

By taking a threshold equal to e = 0.01 we select the spectral components 
corresponding to \c Xn \ ^ for —20 < x n < 140 (see Fig JTB"]) . The error 
between the original and the selected signal is about —18dB. From there the 
spectrum of y(t) splits in three components. 

First component, the reflection on the porthole 

The first component, y~i(t) corresponds to —20 < x n < and is therefore 
defined as: 



m= E <L(v)K(t) 



(34) 



-20 



It is a low frequency signal corresponding to the heterodyne product of the 
probe signal with the reflection on the porthole [22]- It is shown in FigfT41 
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Figure 9: Spectrum of the signal for 9 = tt/5 



Second component, the plasma signal 

The second component has a Fourier spectra that fits the expected be- 
havior corresponding to the reflection of the wave inside the plasma of the 
tokamak, [22]. This component, ^(t), corresponds to < x n < 110 and is 
therefore defined as: 



no 



&(*) = E oxw 



(35) 



X n =0 



It is shown in FigJT5l 

Third component, the multi-reflection component 

The last component corresponds [22] to multi reflections of the waves 
on the wall of the vacuum vessel. This component, ys(t), corresponds to 
110 < x„ < 140 and is therefore defined as: 



no 



ys(t) = £ <L(v)«t) 



(36) 



X n = W 



This component is shown in FigJTBl 

We notice that by undertaking a new factorization of this third component 
it seems possible to separate different successive reflections of the wave but 
this would be out of the scope of this work. 
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Figure 10: Time representation of the reflectometry signal. 



The three components of the reflectometry signal are presented together 
on the same plot (Fig. [T7I) . It is instructive to compare this factorization 
with the original reflectometry signal (see Fig. [TOl) . 



5 Remarks and conclusions 

1) Based on a complete and probabilistically rigorous spectral analysis and 
projection on the eigenvectors of a family of unitary operators, our method 
seems quite robust to disentangle the relevant components of the signals. 
This has been demonstrated both on simulated and on experimental reflec- 
tometry data. In particular in this last case, a clear identification of the 
physical origin of the components and its separation is readily achieved. Such 
separation could not be achieved by the simple filtering techniques. After the 
component separation phase, the method also provides by truncation of some 
subsets of the projection coefficients a very flexible denoising technique. 

Another important conclusion from this study is the fact that by the 
choice of different families of (unitary) operators and their spectral represen- 
tations, different traits and components of the signals may be emphasized. 

2) In the analysis of reflectometry data, component separation and de- 
noising is a required first step to obtain reliable information on the plasma 
density. In particular, accuracy in these measurements is quite critical if in 
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sinus tela 



Figure 11: Tomogram of the reflectometry signal. 



addition to the average local density one also wants to have information on 
plasma fluctuations and turbulence. 

The location x c (u p ) of the reflecting (cut-off) layer for the frequency u p 
is related to the group delay [23] 



d(f) (u>) 
duj 



1 d(j> 



and for a linear frequency sweep of the incident wave 



/(«) = /o + 7* 



one obtains 



df 



1 d<p 
7 ~dt 



(37) 



(38) 



(39) 



Therefore, measurement of the plasma density hinges on an accurate deter- 
mination of the "instantaneous frequency" -$;. Accuracy in the measurement 
of this quantity is quite critical because, the location of the reflecting layer 
being obtained from an integral, errors tend to accumulate. 

Several methods have been devised to obtain the group delay r g from 
the reflectometry data (for a review see [21]). Among them, time-frequency 
analysis [25J was widely explored. The Wigner-Ville (WV) distribution [26] 
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Figure 12: Spectrum c e x of the reflectometry signal y(t) for 9 
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[27], although providing a complete description of the signal in the time- 
frequency plane, raises difficult interpretation problems due to the presence 
of many interference terms that impair the readability of the distribution. 
For this reason the time- frequency method that, so far, has been preferred is 
the spectrogram [23] [2S] |29j . that is, the squared modulus of the short-time 
Fourier transform 



SP (t, f) 



x (u) h{u-t) e- i2nfu du 



h (u) being a peaked short-time window. 

The spectrogram does not really provide the instantaneous frequency, 
because that notion is not well defined anyway All it gives is the product 
of the spectra of x (t) and h (t) . The way the spectrogram is used to infer 
the local rate of phase variation -& is to identify this quantity with the 
maximum or the with the first moment of the spectrogram. An additional 
problem comes about because unwanted phase contributions due to plasma 
turbulence may have a higher amplitude than the contributions due to the 
profile. Correction techniques have been developed to compensate for this 
errors, based for example on Floyd's best path algorithm. The choice of the 
window function is also an important issue and, in particular, an adaptive 
spectrogram technique has been developed to maximize the time-frequency 
concentration [2l] . 

Here, the tomographic method may also provide a needed improvement, 
namely by obtaining the phase derivative directly from the coefficients c n of 
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Figure 13: Part of the spectrum c e Xn {y) of the reflectometry signal used in 
the factorization 

the plasma reflection component. This, together with the selective denoising 
techniques after component separation, will be described in a forthcoming 
paper. 



6 Appendix A. Gauss-Hermite decomposition 
of the tomograms 

From the definition (TTTJT) of the tomogram transform one sees that the cal- 
culation from the data near v = has accuracy problems because of the 
fast variation of the phase in (TlO|) . Two techniques are used to deal with 
this problem. The first one uses a projection of the time signal s (t) on an 
orthogonal basis and the second uses the homogeneity properties (TTTj) and 
an expansion of the Fresnel tomogram near v = 0. The first technique is 
described in this appendix and the second in the Appendix B. 
Let s (t) be a normalized signal 



(40) 



\s(t)\ 2 dt= 1. 
Decompose the signal in Gauss-Hermite polynomials 

oo 

s(t)=J2c n (t)Mt)> 

n=0 
20 



(41) 




Figure 14: First component of the reflectometry signal corresponding to the 
reflection on the porthole 



with 



and 



e -< 2 /2 I 

s(t)ip n (t) dt, 



(42) 
(43) 



Then, the tomogram of the signal is 



M a (X,fi,v) = M Q (X,fi,u) 

with 



E 

n=0 



1 (\ 



T \2 1 -ifi/v 



™/ 2 / ; 

H< • 



2y/k 



M Q {X,n,v) 



V^W+^) 



-xy(^+v*) 



and 



t = t&L. k- n 



IjJL — V 



2 1 — i\ijv 



(44) 
(45) 

(46) 



7 Appendix B. The Fresnel tomogram 

The symplectic tomogram M S (X, //, v) can be reconstructed if one knows the 
(Fresnel) tomogram [30] 
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Figure 15: Second component of the reflectometry signal, corresponding to 
the reflection on the plasma 



M F (X,v) = M,(X,l,u) 

due to the homogeneity property (fTTl) . In fact, one has 



M a (X,n,v) 



1 X v 

— M s -,1,- 



(47) 



(48) 



which means that, if one knows Mf(X,u), the symplectic tomogram is ob- 
tained by replacement and a factor, 



M S (X, //, u) = -^ M F f X -> — ,v -> - 

|/X| \ /i /i 



In terms of signal s(i) it reads: 
M F (A» ' 



2ttM 



D i(X-^) 2 /2^ 



s(y) dy 



'27I-M 



cxp 



— Zf 



i(X-j/) 2 /2^ 



1 d 2 

2 3X 5 



«(j/) ^ 



8{X) 



(49) 



(50) 
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Figure 16: Third component of the reflectometry signal, corresponding to 
multi reflections on the vessel 

Thus for small v one has 



M s (X,/i,z/)^ — 

m 

In the Gauss-Hermite basis it is 






(51) 



M F (X,u) 



,-X 2 /(l+v 2 ) 



v^+ 



V 



ra=0 



1 (I 



\ \2 1-i/u 



nil I 

H< " 



2\fk f 



with 



6=^. i= ^ 



2 — V 

As a series, the Fresnel tomogram is 



2 1-i/u J' 



(52) 
(53) 



M F (X,z/) 



fc=0 



ii/ y i rf 2fe /(^) 



fc! dX 2fe 



(54) 



leading to a symplectic tomogram 



M(X,/i,z/) = ^ 



fc=0 



E £ 



ii/y ld? k f(X/fi,) 



2\i) k\ dX 2k 



(55) 
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Figure 17: The three components of the reflectometry signal. For visual 
purposes, the average of y~i(t) is shifted to 1 and the average of y 3 (t) to -1. 



Since 



1 d 2 t 2 ( 1 

2^ s " (t)+ 2 s " (t) = r + 2 



s n (t), 



n = 0,1,2,.. 



(56) 



on has, for small v, the following Fresnel tomogram of s n (t): 

1 X 2 ' 



M n (X,v = l,v)ns 2 n (X) 



1 



n H i/ 

2 2 



(57) 
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